function [T Y] = problema5(N, infected)

global GAMMA;
global BETA;

BETA = 0.02;
GAMMA = 0.5;

[T, Y] = rk4gral(0, 25, [(1-infected)*N, infected*N, 0], 100, 'fSIR');

plot(T, Y(:,1), 'b');
hold on;
plot(T, Y(:,2), 'r*');
hold on;
plot(T, Y(:,3), 'g+');
hold off;

xlabel('Tiempo en dias (t)');
ylabel('Numero de individuos (N)');
legend(strcat('S(t); S(0) = ',num2str((1-infected)*N)), strcat('I(t); I(0) = ', num2str(infected*N)), 'R(t); R(0) = 0', 'Location', 'NorthEastOutside');
title('Individuos suceptibles, infectados y recuperados en funcion del tiempo en dias');
end
